Bit labeling for amplitude phase shift constellation used with low density parity check (LDPC) codes

ABSTRACT

An approach is provided for bit labeling of a signal constellation. A transmitter generates encoded signals using, according to one embodiment, a structured parity check matrix of a Low Density Parity Check (LDPC) code. The transmitter includes an encoder for transforming an input message into a codeword represented by a plurality of set of bits. The transmitter includes logic for mapping non-sequentially (e.g., interleaving) one set of bits into a higher order constellation (Quadrature Phase Shift Keying (QPSK), 8-PSK, 16-APSK (Amplitude Phase Shift Keying), 32-APSK, etc.), wherein a symbol of the higher order constellation corresponding to the one set of bits is output based on the mapping.

FIELD OF THE INVENTION

[0001] The present invention relates to communication systems, and moreparticularly to coded systems.

BACKGROUND OF THE INVENTION

[0002] Communication systems employ coding to ensure reliablecommunication across noisy communication channels. These communicationchannels exhibit a fixed capacity that can be expressed in terms of bitsper symbol at certain signal to noise ratio (SNR), defining atheoretical upper limit (known as the Shannon limit). As a result,coding design has aimed to achieve rates approaching this Shannon limit.Conventional coded communication systems have separately treated theprocesses of coding and modulation. Moreover, little attention has beenpaid to labeling of signal constellations.

[0003] A signal constellation provides a set of possible symbols thatare to be transmitted, whereby the symbols correspond to codewordsoutput from an encoder. One choice of constellation labeling involvesGray-code labeling. With Gray-code labeling, neighboring signal pointsdiffer in exactly one bit position. The prevailing conventional view ofmodulation dictates that any reasonable labeling scheme can be utilized,which in part is responsible for the paucity of research in this area.

[0004] With respect to coding, one class of codes that approach theShannon limit is Low Density Parity Check (LDPC) codes. Traditionally,LDPC codes have not been widely deployed because of a number ofdrawbacks. One drawback is that the LDPC encoding technique is highlycomplex. Encoding an LDPC code using its generator matrix would requirestoring a very large, non-sparse matrix. Additionally, LDPC codesrequire large blocks to be effective; consequently, even though paritycheck matrices of LDPC codes are sparse, storing these matrices isproblematic.

[0005] From an implementation perspective, a number of challenges areconfronted. For example, storage is an important reason why LDPC codeshave not become widespread in practice. Also, a key challenge in LDPCcode implementation has been how to achieve the connection networkbetween several processing engines (nodes) in the decoder. Further, thecomputational load in the decoding process, specifically the check nodeoperations, poses a problem.

[0006] Therefore, there is a need for a bit labeling approach thatsupplements code performance of coded systems in general. There is alsoa need for using LDPC codes efficiently to support high data rates,without introducing greater complexity. There is also a need to improveperformance of LDPC encoders and decoders.

SUMMARY OF THE INVENTION

[0007] These and other needs are addressed by the present invention,wherein an approach is provided for bit labeling of a signalconstellation. An encoder, such as a Low Density Parity Check (LDPC)encoder, generates encoded signals by transforming an input message intoa codeword represented by a plurality of set of bits. These bits aremapped non-sequentially (e.g., interleaving) a higher orderconstellation (Quadrature Phase Shift Keying (QPSK), 8-PSK, 16-APSK(Amplitude Phase Shift Keying), 32-APSK, etc. The above arrangementadvantageously provides enhanced performance of the codes.

[0008] According to one aspect of an embodiment of the presentinvention, a method for transmitting encoded signals is disclosed. Themethod includes receiving one of a plurality of set of bits of acodeword from an encoder for transforming an input message into thecodeword. The method also includes non-sequentially mapping the one setof bits into a higher order constellation. Further, the method includesoutputting a symbol of the higher order constellation corresponding tothe one set of bits based on the mapping.

[0009] According to another aspect of an embodiment of the presentinvention, a transmitter for generating encoded signals is disclosed.The transmitter includes an encoder configured to transform an inputmessage into a codeword represented by a plurality of set of bits.Additionally, the transmitter includes logic configured to mapnon-sequentially one set of bits into a higher order constellation,wherein a symbol of the higher order constellation corresponding to theone set of bits is output based on the mapping.

[0010] According to another aspect of an embodiment of the presentinvention, a method for processing encoded signals is disclosed. Themethod includes demodulating a received encoded signal representing acodeword, wherein the encoded signal being modulated according to anon-sequential mapping of a plurality of bits corresponding to thecodeword. The method also includes decoding the codeword associated withthe encoded signal.

[0011] Still other aspects, features, and advantages of the presentinvention are readily apparent from the following detailed description,simply by illustrating a number of particular embodiments andimplementations, including the best mode contemplated for carrying outthe present invention. The present invention is also capable of otherand different embodiments, and its several details can be modified invarious obvious respects, all without departing from the spirit andscope of the present invention. Accordingly, the drawing and descriptionare to be regarded as illustrative in nature, and not as restrictive.

BRIEF DESCRIPTION OF THE DRAWINGS

[0012] The present invention is illustrated by way of example, and notby way of limitation, in the figures of the accompanying drawings and inwhich like reference numerals refer to similar elements and in which:

[0013]FIG. 1 is a diagram of a communications system configured toutilize Low Density Parity Check (LDPC) codes, according to anembodiment of the present invention;

[0014]FIGS. 2A and 2B are diagrams of exemplary LDPC encoders deployedin the transmitter of FIG. 1;

[0015]FIG. 3 is a diagram of an exemplary receiver in the system ofFIG., 1;

[0016]FIG. 4 is a diagram of a sparse parity check matrix, in accordancewith an embodiment of the present invention;

[0017]FIG. 5 is a diagram of a bipartite graph of an LDPC code of thematrix of FIG. 4;

[0018]FIG. 6 is a diagram of a sub-matrix of a sparse parity checkmatrix, wherein the sub-matrix contains parity check values restrictedto the lower triangular region, according to an embodiment of thepresent invention;

[0019]FIG. 7 is a graph showing performance between codes utilizingunrestricted parity check matrix (H matrix) versus restricted H matrixhaving a sub-matrix as in FIG. 6;

[0020]FIGS. 8A and 8B are, respectively, a diagram of a non-Gray 8-PSKmodulation scheme, and a Gray 8-PSK modulation, each of which can beused in the system of FIG. 1;

[0021]FIG. 8C is a diagram of a process for bit labeling for a higherorder signal constellation, according to an embodiment of the presentinvention;

[0022]FIG. 8D is a diagram of exemplary 16-APSK (Amplitude Phase ShiftKeying) constellations;

[0023]FIG. 8E is a graph of Packet Error Rate (PER) versussignal-to-noise for the constellations of FIG. 8D;

[0024]FIG. 8F is a diagram of constellations for Quadrature Phase ShiftKeying (QPSK), 8-PSK, 16-APSK and 32-APSK symbols, in accordance with anembodiment of the present invention;

[0025]FIG. 8G is a diagram of alternative constellations for 8-PSK,16-APSK and 32-APSK symbols, in accordance with an embodiment of thepresent invention;

[0026]FIG. 8H is a graph of Packet Error Rate (PER) versussignal-to-noise for the constellations of FIG. 8F;

[0027]FIG. 9 is a graph showing performance between codes utilizing Graylabeling versus non-Gray labeling;

[0028]FIG. 10 is a flow chart of the operation of the LDPC decoder usingnon-Gray mapping, according to an embodiment of the present invention;

[0029]FIG. 11 is a flow chart of the operation of the LDPC decoder ofFIG. 3 using Gray mapping, according to an embodiment of the presentinvention;

[0030] FIGS. 12A-12C are diagrams of the interactions between the checknodes and the bit nodes in a decoding process, according to anembodiment of the present invention;

[0031]FIGS. 13A and 13B are flowcharts of processes for computingoutgoing messages between the check nodes and the bit nodes using,respectively, a forward-backward approach and a parallel approach,according to various embodiments of the present invention;

[0032] FIGS. 14A-14 are graphs showing simulation results of LDPC codesgenerated in accordance with various embodiments of the presentinvention;

[0033]FIGS. 15A and 15B are diagrams of the top edge and bottom edge,respectively, of memory organized to support structured access as torealize randomness in LDPC coding, according to an embodiment of thepresent invention; and

[0034]FIG. 16 is a diagram of a computer system that can perform theprocesses of encoding and decoding of LDPC codes, in accordance withembodiments of the present invention.

DESCRIPTION OF THE PREFERRED EMBODIMENT

[0035] A system, method, and software for bit labeling for signalconstellations are described. In the following description, for thepurposes of explanation, numerous specific details are set forth inorder to provide a thorough understanding of the present invention. Itis apparent, however, to one skilled in the art that the presentinvention may be practiced without these specific details or with anequivalent arrangement. In other instances, well-known structures anddevices are shown in block diagram form in order to avoid unnecessarilyobscuring the present invention.

[0036] Although the present invention is described with respect to LDPCcodes, it is recognized that the bit labeling approach can be utilizedwith other codes. Further, this approach can be implemented with uncodedsystems.

[0037]FIG. 1 is a diagram of a communications system configured toutilize Low Density Parity Check (LDPC) codes, according to anembodiment of the present invention. A digital communications system 100includes a transmitter 101 that generates signal waveforms across acommunication channel 103 to a receiver 105. In this discretecommunications system 100, the transmitter 101 has a message source thatproduces a discrete set of possible messages; each of the possiblemessages has a corresponding signal waveform. These signal waveforms areattenuated, or otherwise altered, by communications channel 103. Tocombat the noise channel 103, LDPC codes are utilized.

[0038] The LDPC codes that are generated by the transmitter 101 enableshigh speed implementation without incurring any performance loss. Thesestructured LDPC codes output from the transmitter 101 avoid assignmentof a small number of check nodes to the bit nodes already vulnerable tochannel errors by virtue of the modulation scheme (e.g., 8-PSK).

[0039] Such LDPC codes have a parallelizable decoding algorithm (unliketurbo codes), which advantageously involves simple operations such asaddition, comparison and table look-up. Moreover, carefully designedLDPC codes do not exhibit any sign of error floor.

[0040] According to one embodiment of the present invention, thetransmitter 101 generates, using a relatively simple encoding technique,LDPC codes based on parity check matrices (which facilitate efficientmemory access during decoding) to communicate with the receiver 105. Thetransmitter 101 employs LDPC codes that can outperform concatenatedturbo+RS (Reed-Solomon) codes, provided the block length is sufficientlylarge.

[0041]FIGS. 2A and 2B are diagrams of exemplary LDPC encoders deployedin the transmitter of FIG. 1. As seen in FIG. 2A, a transmitter 200 isequipped with an LDPC encoder 203 that accepts input from an informationsource 201 and outputs coded stream of higher redundancy suitable forerror correction processing at the receiver 105. The information source201 generates k signals from a discrete alphabet, X. LDPC codes arespecified with parity check matrices. On the other hand, encoding LDPCcodes require, in general, specifying the generator matrices. Eventhough it is possible to obtain generator matrices from parity checkmatrices using Gaussian elimination, the resulting matrix is no longersparse and storing a large generator matrix can be complex.

[0042] Encoder 203 generates signals from alphabet Y to a signal mapper206, which provides a mapping of the alphabet Y to the symbols of thesignal constellation corresponding to the modulation scheme employed bya modulator 205. This mapping, according to one embodiment of thepresent invention, follows a non-sequential scheme, such asinterleaving. Exemplary mappings are more fully described below withrespect to FIG. 8C. The encoder 203 uses a simple encoding techniquethat makes use of only the parity check matrix by imposing structureonto the parity check matrix. Specifically, a restriction is placed onthe parity check matrix by constraining certain portion of the matrix tobe triangular. The construction of such a parity check matrix isdescribed more fully below in FIG. 6. Such a restriction results innegligible performance loss, and therefore, constitutes an attractivetradeoff.

[0043] The modulator 205 modulates the symbols of the signalconstellation from the mapper 206 to signal waveforms that aretransmitted to a transmit antenna 207, which emits these waveforms overthe communication channel 103. The transmissions from the transmitantenna 207 propagate to a receiver, as discussed below.

[0044]FIG. 2B shows an LDPC encoder utilized with a Bose ChaudhuriHocquenghem (BCH) encoder and a cyclic redundancy check (CRC) encoder,according to one embodiment of the present invention. Under thisscenario, the codes generated by the LDPC encoder 203, along with theCRC encoder 209 and the BCH encoder 211, have a concatenated outer BCHcode and inner low density parity check (LDPC) code. Furthermore, errordetection is achieved using cyclic redundancy check (CRC) codes. The CRCencoder 209, in an exemplary embodiment, encodes using an 8-bit CRC codewith generator polynomial (x⁵+x⁴+x³+x²+1)(x²+x+1)(x+1).

[0045] The LDPC encoder 203 systematically encodes an information blockof size k_(ldpc), i=(i₀, i₁, . . . , i_(k) _(ldpc) ⁻¹) onto a codewordof size n_(ldpc), c=(i₀, i₁, . . . , i_(k) _(ldpc) ⁻¹, p₀, p₁, . . .p_(n) _(ldpc) _(−k) _(ldpc) ⁻¹) The transmission of the codeword startsin the given order from i₀ and ends with p_(n) _(ldpc) _(−k) _(ldpc) ⁻¹.LDPC code parameters (n_(ldpc), k_(ldpc)) are given in Table 1 below.TABLE 1 LDPC Code Parameters (n_(ldpc), k_(ldpc)) LDPC Uncoded LDPCCoded Block Block Length Length Code Rate k_(ldpc) n_(ldpc) 1/2 3240064800 2/3 43200 64800 3/4 48600 64800 4/5 51840 64800 5/6 54000 648003/5 38880 64800 8/9 57600 64800  9/10 58320 64800

[0046] The task of the LDPC encoder 203 is to determinen_(ldpc)−k_(ldpc) parity bits (p₀, p₁, . . . , p_(n) _(ldpc) _(−k)_(ldpc) ⁻¹) for every block of k_(ldpc) information bits, (i₀, i₁, . . ., i_(k) _(ldpc) ⁻¹). The procedure is as follows. First, the parity bitsare initialized; p₀=p₁=p₂= . . . =p_(n) _(ldpc) _(−k) _(ldpc) ⁻¹=0. Thefirst information bit, i₀, are accumulated at parity bit addressesspecified in the first row of Tables 3 through 10. For example, for rate⅔ (Table 3), the following results:

[0047] p₀=p₀⊕i₀

[0048] p₁₀₄₉₁=p₁₀₄₉₁⊕i₀

[0049] p₁₆₀₄₃=p₁₆₀₄₃⊕i₀

[0050] p₅₀₆=p₅₀₆⊕i₀

[0051] p₁₂₈₂₆=p₁₂₈₂₆⊕i₀

[0052] p₈₀₆₅=p₈₀₆₅⊕i₀

[0053] p₈₂₂₆=p₈₂₂₆⊕i₀

[0054] p₂₇₆₇=p₂₇₆₇⊕i₀

[0055] p₂₄₀=p₂₄₀⊕i₀

[0056] p₁₈₆₇₃=p₁₈₆₇₃⊕i₀

[0057] p₉₂₇₉=p₉₂₇₉⊕i₀

[0058] p₁₀₅₇₉=p₁₀₅₇₉⊕i₀

[0059] p₂₀₉₂₈=p₂₀₉₂₈⊕i₀

[0060] (All additions are in GF(2)).

[0061] Then, for the next 359 information bits, i_(m), m=1, 2, . . . ,359, these bits are accumulated at parity bit addresses {x+m mod360×q}mod(n_(ldpc)−k_(ldpc)), where x denotes the address of the paritybit accumulator corresponding to the first bit i₀, and q is a code ratedependent constant specified in Table 2. Continuing with the example,q=60 for rate ⅔. By way of example, for information bit i₁, thefollowing operations are performed:

[0062] p₆₀=p₆₀⊕i₁

[0063] p₁₀₅₅₁=p₁₀₅₅₁⊕i₁

[0064] p₁₆₁₀₃=p₁₆₁₀₃⊕i₁

[0065] p₅₆₆=p₅₆₆⊕i₁

[0066] p₁₂₈₈₆=p₁₂₈₈₆⊕i₁

[0067] p₈₁₂₅=p₈₁₂₅⊕i₁

[0068] p₈₂₈₆=p₈₂₈₆⊕i₁

[0069] p₂₈₂₇=p₂₈₂₇⊕i₁

[0070] p₃₀₀=p₃₀₀⊕i₁

[0071] p₁₈₇₃₃=p₁₈₇₃₃⊕i₁

[0072] p₉₃₃₉=p₉₃₃₉⊕i₁

[0073] p₁₀₆₃₉=p₁₀₆₃₉⊕i₁

[0074] p₂₀₉₈₈=p₂₀₉₈₈⊕i₁

[0075] For the 361^(st) information bit i₃₆₀, the addresses of theparity bit accumulators are given in the second row of the Tables 3through 10. In a similar manner the addresses of the parity bitaccumulators for the following 359 information bits i_(m), m=361, 362, .. . , 719 are obtained using the formula {x+m mod 360×q}mod(n_(ldpc)−k_(ldpc)), where x denotes the address of the parity bitaccumulator corresponding to the information bit i₃₆₀, i.e., the entriesin the second row of the Tables 3-10. In a similar manner, for everygroup of 360 new information bits, a new row from Tables 3 through 10are used to find the addresses of the parity bit accumulators.

[0076] After all of the information bits are exhausted, the final paritybits are obtained as follows. First, the following operations areperformed, starting with i=1

p _(i)=p_(i)⊕p_(i−1), i=1, 2, . . . , n_(ldpc)−k_(ldpc)−1.

[0077] Final content of p_(i), i=0,1, . . . , n_(ldpc)−k_(ldpc)−1 isequal to the parity bit p_(i). TABLE 2 Code Rate q 2/3 60 5/6 30 1/2 903/4 45 4/5 36 3/5 72 8/9 20  9/10 18

[0078] TABLE 3 Address of Parity Bit Accumulators (Rate 2/3) 0 1049116043 506 12826 8065 8226 2767 240 18673 9279 10579 20928 1 17819 83136433 6224 5120 5824 12812 17187 9940 13447 13825 18483 2 17957 6024 868118628 12794 5915 14576 10970 12064 20437 4455 7151 3 19777 6183 997214536 8182 17749 11341 5556 4379 17434 15477 18532 4 4651 19689 1608 65916707 14335 6143 3058 14618 17894 20684 5306 5 9778 2552 12096 1236915198 16890 4851 3109 1700 18725 1997 15882 6 486 6111 13743 11537 55917433 15227 14145 1483 3887 17431 12430 7 20647 14311 11734 4180 81105525 12141 15761 18661 18441 10569 8192 8 3791 14759 15264 19918 101329062 10010 12786 10675 9682 19246 5454 9 19525 9485 7777 19999 8378 92093163 20232 6690 16518 716 7353 10 4588 6709 20202 10905 915 4317 1107313576 16433 368 3508 21171 11 14072 4033 19959 12608 631 19494 141608249 10223 21504 12395 4322 12 13800 14161 13 2948 9647 14 14693 1602715 20506 11082 16 1143 9020 17 13501 4014 18 1548 2190 19 12216 21556 202095 19897 21 4189 7958 22 15940 10048 23 515 12614 24 8501 8450 2517595 16784 26 5913 8495 27 16394 10423 28 7409 6981 29 6678 15939 3020344 12987 31 2510 14588 32 17918 6655 33 6703 19451 34 496 4217 357290 5766 36 10521 8925 37 20379 11905 38 4090 5838 39 19082 17040 4020233 12352 41 19365 19546 42 6249 19030 43 11037 19193 44 19760 1177245 19644 7428 46 16076 3521 47 11779 21062 48 13062 9682 49 8934 5217 5011087 3319 51 18892 4356 52 7894 3898 53 5963 4360 54 7346 11726 55 51825609 56 2412 17295 57 9845 20494 58 6687 1864 59 20564 5216 0 1822617207 1 9380 8266 2 7073 3065 3 18252 13437 4 9161 15642 5 10714 10153 611585 9078 7 5359 9418 8 9024 9515 9 1206 16354 10 14994 1102 11 937520796 12 15964 6027 13 14789 6452 14 8002 18591 15 14742 14089 16 2533045 17 1274 19286 18 14777 2044 19 13920 9900 20 452 7374 21 18206 992122 6131 5414 23 10077 9726 24 12045 5479 25 4322 7990 26 15616 5550 2715561 10661 28 20718 7387 29 2518 18804 30 8984 2600 31 6516 17909 3211148 98 33 20559 3704 34 7510 1569 35 16000 11692 36 9147 10303 3716650 191 38 15577 18685 39 17167 20917 40 4256 3391 41 20092 17219 429218 5056 43 18429 8472 44 12093 20753 45 16345 12748 46 16023 11095 475048 17595 48 18995 4817 49 16483 3536 50 1439 16148 51 3661 3039 5219010 18121 53 8968 11793 54 13427 18003 55 5303 3083 56 531 16668 574771 6722 58 5695 7960 59 3589 14630

[0079] TABLE 4 Address of Parity Bit Accumulators (Rate 5/6) 0 4362 4168909 4156 3216 3112 2560 2912 6405 8593 4969 6723 1 2479 1786 8978 30114339 9313 6397 2957 7288 5484 6031 10217 2 10175 9009 9889 3091 49857267 4092 8874 5671 2777 2189 8716 3 9052 4795 3924 3370 10058 1128 999610165 9360 4297 434 5138 4 2379 7834 4835 2327 9843 804 329 8353 71673070 1528 7311 5 3435 7871 348 3693 1876 6585 10340 7144 5870 2084 40522780 6 3917 3111 3476 1304 10331 5939 5199 1611 1991 699 8316 9960 76883 3237 1717 10752 7891 9764 4745 3888 10009 4176 4614 1567 8 105872195 1689 2968 5420 2580 2883 6496 111 6023 1024 4449 9 3786 8593 20743321 5057 1450 3840 5444 6572 3094 9892 1512 10 8548 1848 10372 45857313 6536 6379 1766 9462 2456 5606 9975 11 8204 10593 7935 3636 3882 3945968 8561 2395 7289 9267 9978 12 7795 74 1633 9542 6867 7352 6417 756810623 725 2531 9115 13 7151 2482 4260 5003 10105 7419 9203 6691 87982092 8263 3755 14 3600 570 4527 200 9718 6771 1995 8902 5446 768 11036520 15 6304 7621 16 6498 9209 17 7293 6786 18 5950 1708 19 8521 1793 206174 7854 21 9773 1190 22 9517 10268 23 2181 9349 24 1949 5560 25 1556555 26 8600 3827 27 5072 1057 28 7928 3542 29 3226 3762 0 7045 2420 19645 2641 2 2774 2452 3 5331 2031 4 9400 7503 5 1850 2338 6 10456 9774 71692 9276 8 10037 4038 9 3964 338 10 2640 5087 11 858 3473 12 5582 568313 9523 916 14 4107 1559 15 4506 3491 16 8191 4182 17 10192 6157 18 56683305 19 3449 1540 20 4766 2697 21 4069 6675 22 1117 1016 23 5619 3085 248483 8400 25 8255 394 26 6338 5042 27 6174 5119 28 7203 1989 29 17815174 0 1464 3559 1 3376 4214 2 7238 67 3 10595 8831 4 1221 6513 5 53004652 6 1429 9749 7 7878 5131 8 4435 10284 9 6331 5507 10 6662 4941 119614 10238 12 8400 8025 13 9156 5630 14 7067 8878 15 9027 3415 16 16903866 17 2854 8469 18 6206 630 19 363 5453 20 4125 7008 21 1612 6702 229069 9226 23 5767 4060 24 3743 9237 25 7018 5572 26 8892 4536 27 8536064 28 8069 5893 29 2051 2885 0 10691 3153 1 3602 4055 2 328 1717 32219 9299 4 1939 7898 5 617 206 6 8544 1374 7 10676 3240 8 6672 9489 93170 7457 10 7868 5731 11 6121 10732 12 4843 9132 13 580 9591 14 62679290 15 3009 2268 16 195 2419 17 8016 1557 18 1516 9195 19 8062 9064 202095 8968 21 753 7326 22 6291 3833 23 2614 7844 24 2303 646 25 2075 61126 4687 362 27 8684 9940 28 4830 2065 29 7038 1363 0 1769 7837 1 38011689 2 10070 2359 3 3667 9918 4 1914 6920 5 4244 5669 6 10245 7821 77648 3944 8 3310 5488 9 6346 9666 10 7088 6122 11 1291 7827 12 105928945 13 3609 7120 14 9168 9112 15 6203 8052 16 3330 2895 17 4264 1056318 10556 6496 19 8807 7645 20 1999 4530 21 9202 6818 22 3403 1734 232106 9023 24 6881 3883 25 3895 2171 26 4062 6424 27 3755 9536 28 46832131 29 7347 8027

[0080] TABLE 5 Address of Parity Bit Accumulators (Rate 1/2) 54 931814392 27561 26909 10219 2534 8597 55 7263 4635 2530 28130 3033 238303651 56 24731 23583 26036 17299 5750 792 9169 57 5811 26154 18653 1155115447 13685 16264 58 12610 11347 28768 2792 3174 29371 12997 59 1678916018 21449 6165 21202 15850 3186 60 31016 21449 17618 6213 12166 833418212 61 22836 14213 11327 5896 718 11727 9308 62 2091 24941 29966 236349013 15587 5444 63 22207 3983 16904 28534 21415 27524 25912 64 256874501 22193 14665 14798 16158 5491 65 4520 17094 23397 4264 22370 1694121526 66 10490 6182 32370 9597 30841 25954 2762 67 22120 22865 2987015147 13668 14955 19235 68 6689 18408 18346 9918 25746 5443 20645 6929982 12529 13858 4746 30370 10023 24828 70 1262 28032 29888 13063 2403321951 7863 71 6594 29642 31451 14831 9509 9335 31552 72 1358 6454 1663320354 24598 624 5265 73 19529 295 18011 3080 13364 8032 15323 74 119811510 7960 21462 9129 11370 25741 75 9276 29656 4543 30699 20646 2192128050 76 15975 25634 5520 31119 13715 21949 19605 77 18688 4608 3175530165 13103 10706 29224 78 21514 23117 12245 26035 31656 25631 30699 799674 24966 31285 29908 17042 24588 31857 80 21856 27777 29919 2700014897 11409 7122 81 29773 23310 263 4877 28622 20545 22092 82 15605 565121864 3967 14419 22757 15896 83 30145 1759 10139 29223 26086 10556 509884 18815 16575 2936 24457 26738 6030 505 85 30326 22298 27562 2013126390 6247 24791 86 928 29246 21246 12400 15311 32309 18608 87 203146025 26689 16302 2296 3244 19613 88 6237 11943 22851 15642 23857 1511220947 89 26403 25168 19038 18384 8882 12719 7093 0 14567 24965 1 3908100 2 10279 240 3 24102 764 4 12383 4173 5 13861 15918 6 21327 1046 75288 14579 8 28158 8069 9 16583 11098 10 16681 28363 11 13980 24725 1232169 17989 13 10907 2767 14 21557 3818 15 26676 12422 16 7676 8754 1714905 20232 18 15719 24646 19 31942 8589 20 19978 27197 21 27060 1507122 6071 26649 23 10393 11176 24 9597 13370 25 7081 17677 26 1433 1951327 26925 9014 28 19202 8900 29 18152 30647 30 20803 1737 31 11804 2522132 31683 17783 33 29694 9345 34 12280 26611 35 6526 26122 36 26165 1124137 7666 26962 38 16290 8480 39 11774 10120 40 30051 30426 41 1335 1542442 6865 17742 43 31779 12489 44 32120 21001 45 14508 6996 46 979 2502447 4554 21896 48 7989 21777 49 4972 20661 50 6612 2730 51 12742 4418 5229194 595 53 19267 20113

[0081] TABLE 6 Address of Parity Bit Accumulators (Rate 3/4) 0 6385 790114611 13389 11200 3252 5243 2504 2722 821 7374 1 11359 2698 357 1382412772 7244 6752 15310 852 2001 11417 2 7862 7977 6321 13612 12197 1444915137 13860 1708 6399 13444 3 1560 11804 6975 13292 3646 3812 8772 73065795 14327 7866 4 7626 11407 14599 9689 1628 2113 10809 9283 1230 152414870 5 1610 5699 15876 9446 12515 1400 6303 5411 14181 13925 7358 6 40598836 3405 7853 7992 15336 5970 10368 10278 9675 4651 7 4441 3963 91532109 12683 7459 12030 12221 629 15212 406 8 6007 8411 5771 3497 54314202 875 9186 6235 13908 3563 9 3232 6625 4795 546 9781 2071 7312 33997250 4932 12652 10 8820 10088 11090 7069 6585 13134 10158 7183 488 74559238 11 1903 10818 119 215 7558 11046 10615 11545 14784 7961 15619 123655 8736 4917 15874 5129 2134 15944 14768 7150 2692 1469 13 8316 3820505 8923 6757 806 7957 4216 15589 13244 2622 14 14463 4852 15733 304111193 12860 13673 8152 6551 15108 8758 15 3149 11981 16 13416 6906 1713098 13352 18 2009 14460 19 7207 4314 20 3312 3945 21 4418 6248 22 266913975 23 7571 9023 24 14172 2967 25 7271 7138 26 6135 13670 27 749014559 28 8657 2466 29 8599 12834 30 3470 3152 31 13917 4365 32 602413730 33 10973 14182 34 2464 13167 35 5281 15049 36 1103 1849 37 20581069 38 9654 6095 39 14311 7667 40 15617 8146 41 4588 11218 42 136606243 43 8578 7874 44 11741 2686 0 1022 1264 1 12604 9965 2 8217 2707 33156 11793 4 354 1514 5 6978 14058 6 7922 16079 7 15087 12138 8 50536470 9 12687 14932 10 15458 1763 11 8121 1721 12 12431 549 13 4129 709114 1426 8415 15 9783 7604 16 6295 11329 17 1409 12061 18 8065 9087 192918 8438 20 1293 14115 21 3922 13851 22 3851 4000 23 5865 1768 24 265514957 25 5565 6332 26 4303 12631 27 11653 12236 28 16025 7632 29 465514128 30 9584 13123 31 13987 9597 32 15409 12110 33 8754 15490 34 741615325 35 2909 15549 36 2995 8257 37 9406 4791 38 11111 4854 39 2812 852140 8476 14717 41 7820 15360 42 1179 7939 43 2357 8678 44 7703 6216 03477 7067 1 3931 13845 2 7675 12899 3 1754 8187 4 7785 1400 5 9213 58916 2494 7703 7 2576 7902 8 4821 15682 9 10426 11935 10 1810 904 11 113329264 12 11312 3570 13 14916 2650 14 7679 7842 15 6089 13084 16 3938 275117 8509 4648 18 12204 8917 19 5749 12443 20 12613 4431 21 1344 4014 228488 13850 23 1730 14896 24 14942 7126 25 14983 8863 26 6578 8564 274947 396 28 297 12805 29 13878 6692 30 11857 11186 31 14395 11493 3216145 12251 33 13462 7428 34 14526 13119 35 2535 11243 36 6465 12690 376872 9334 38 15371 14023 39 8101 10187 40 11963 4848 41 15125 6119 428051 14465 43 11139 5167 44 2883 14521

[0082] TABLE 7 Address of Parity Bit Accumulators (Rate 4/5) 0 149 112125575 6360 12559 8108 8505 408 10026 12828 1 5237 490 10677 4998 38693734 3092 3509 7703 10305 2 8742 5553 2820 7085 12116 10485 564 77952972 2157 3 2699 4304 8350 712 2841 3250 4731 10105 517 7516 4 120671351 11992 12191 11267 5161 537 6166 4246 2363 5 6828 7107 2127 37245743 11040 10756 4073 1011 3422 6 11259 1216 9526 1466 10816 940 37442815 11506 11573 7 4549 11507 1118 1274 11751 5207 7854 12803 4047 64848 8430 4115 9440 413 4455 2262 7915 12402 8579 7052 9 3885 9126 56654505 2343 253 4707 3742 4166 1556 10 1704 8936 6775 8639 8179 7954 82347850 8883 8713 11 11716 4344 9087 11264 2274 8832 9147 11930 6054 545512 7323 3970 10329 2170 8262 3854 2087 12899 9497 11700 13 4418 14672490 5841 817 11453 533 11217 11962 5251 14 1541 4525 7976 3457 95367725 3788 2982 6307 5997 15 11484 2739 4023 12107 6516 551 2572 66288150 9852 16 6070 1761 4627 6534 7913 3730 11866 1813 12306 8249 1712441 5489 8748 7837 7660 2102 11341 2936 6712 11977 18 10155 4210 191010 10483 20 8900 10250 21 10243 12278 22 7070 4397 23 12271 3887 2411980 6836 25 9514 4356 26 7137 10281 27 11881 2526 28 1969 11477 293044 10921 30 2236 8724 31 9104 6340 32 7342 8582 33 11675 10405 34 646712775 35 3186 12198 0 9621 11445 1 7486 5611 2 4319 4879 3 2196 344 47527 6650 5 10693 2440 6 6755 2706 7 5144 5998 8 11043 8033 9 4846 443510 4157 9228 11 12270 6562 12 11954 7592 13 7420 2592 14 8810 9636 15689 5430 16 920 1304 17 1253 11934 18 9559 6016 19 312 7589 20 4439 419721 4002 9555 22 12232 7779 23 1494 8782 24 10749 3969 25 4368 3479 266316 5342 27 2455 3493 28 12157 7405 29 6598 11495 30 11805 4455 31 96252090 32 4731 2321 33 3578 2608 34 8504 1849 35 4027 1151 0 5647 4935 14219 1870 2 10968 8054 3 6970 5447 4 3217 5638 5 8972 669 6 5618 12472 71457 1280 8 8868 3883 9 8866 1224 10 8371 5972 11 266 4405 12 3706 324413 6039 5844 14 7200 3283 15 1502 11282 16 12318 2202 17 4523 965 189587 7011 19 2552 2051 20 12045 10306 21 11070 5104 22 6627 6906 23 98892121 24 829 9701 25 2201 1819 26 6689 12925 27 2139 8757 28 12004 594829 8704 3191 30 8171 10933 31 6297 7116 32 616 7146 33 5142 9761 3410377 8138 35 7616 5811 0 7285 9863 1 7764 10867 2 12343 9019 3 44148331 4 3464 642 5 6960 2039 6 786 3021 7 710 2086 8 7423 5601 9 81204885 10 12385 11990 11 9739 10034 12 424 10162 13 1347 7597 14 1450 11215 7965 8478 16 8945 7397 17 6590 8316 18 6838 9011 19 6174 9410 20 255113 21 6197 5835 22 12902 3844 23 4377 3505 24 5478 8672 25 4453 2132 269724 1380 27 12131 11526 28 12323 9511 29 8231 1752 30 497 9022 31 92883080 32 2481 7515 33 2696 268 34 4023 12341 35 7108 5553

[0083] TABLE 8 Address of Parity Bit Accumulators (Rate 3/5) 22422 1028211626 19997 11161 2922 3122 99 5625 17064 8270 179 25087 16218 17015 82820041 25656 4186 11629 22599 17305 22515 6463 11049 22853 25706 143885500 19245 8732 2177 13555 11346 17265 3069 16581 22225 12563 1971723577 11555 25496 6853 25403 5218 15925 21766 16529 14487 7643 1071517442 11119 5679 14155 24213 21000 1116 15620 5340 8636 16693 1434 56356516 9482 20189 1066 15013 25361 14243 18506 22236 20912 8952 5421 156916126 21595 500 6904 13059 6802 8433 4694 5524 14216 3685 19721 254209937 23813 9047 25651 16826 21500 24814 6344 17382 7064 13929 4004 1655212818 8720 5286 2206 22517 2429 19065 2921 21611 1873 7507 5661 2300623128 20543 19777 1770 4636 20900 14931 9247 12340 11008 12966 4471 273116445 791 6635 14556 18865 22421 22124 12697 9803 25485 7744 18254 113139004 19982 23963 18912 7206 12500 4382 20067 6177 21007 1195 23547 24837756 11158 14646 20534 3647 17728 11676 11843 12937 4402 8261 22944 930624009 10012 11081 3746 24325 8060 19826 842 8836 2898 5019 7575 745525244 4736 14400 22981 5543 8006 24203 13053 1120 5128 3482 9270 1305915825 7453 23747 3656 24585 16542 17507 22462 14670 15627 15290 419822748 5842 13395 23918 16985 14929 3726 25350 24157 24896 16365 1642313461 16615 8107 24741 3604 25904 8716 9604 20365 3729 17245 18448 986220831 25326 20517 24618 13282 5099 14183 8804 16455 17646 15376 1819425528 1777 6066 21855 14372 12517 4488 17490 1400 8135 23375 20879 84764084 12936 25536 22309 16582 6402 24360 25119 23586 128 4761 10443 225368607 9752 25446 15053 1856 4040 377 21160 13474 5451 17170 5938 1025611972 24210 17833 22047 16108 13075 9648 24546 13150 23867 7309 197982988 16858 4825 23950 15125 20526 3553 11525 23366 2452 17626 1926520172 18060 24593 13255 1552 18839 21132 20119 15214 14705 7096 101745663 18651 19700 12524 14033 4127 2971 17499 16287 22368 21463 794318880 5567 8047 23363 6797 10651 24471 14325 4081 7258 4949 7044 1078797 22910 20474 4318 21374 13231 22985 5056 3821 23718 14178 9978 1903023594 8895 25358 6199 22056 7749 13310 3999 23697 16445 22636 5225 2243724153 9442 7978 12177 2893 20778 3175 8645 11863 24623 10311 25767 170573691 20473 11294 9914 22815 2574 8439 3699 5431 24840 21908 16088 182448208 5755 19059 8541 24924 6454 11234 10492 16406 10831 11436 9649 1626411275 24953 2347 12667 19190 7257 7174 24819 2938 2522 11749 3627 596913862 1538 23176 6353 2855 17720 2472 7428 573 15036 0 18539 18661 110502 3002 2 9368 10761 3 12299 7828 4 15048 13362 5 18444 24640 6 2077519175 7 18970 10971 8 5329 19982 9 11296 18655 10 15046 20659 11 730022140 12 22029 14477 13 11129 742 14 13254 13813 15 19234 13273 16 607921122 17 22782 5828 18 19775 4247 19 1660 19413 20 4403 3649 21 1337125851 22 22770 21784 23 10757 14131 24 16071 21617 25 6393 3725 26 59719968 27 5743 8084 28 6770 9548 29 4285 17542 30 13568 22599 31 17864617 32 23238 11648 33 19627 2030 34 13601 13458 35 13740 17328 36 2501213944 37 22513 6687 38 4934 12587 39 21197 5133 40 22705 6938 41 753424633 42 24400 12797 43 21911 25712 44 12039 1140 45 24306 1021 46 1401220747 47 11265 15219 48 4670 15531 49 9417 14359 50 2415 6504 51 2496424690 52 14443 8816 53 6926 1291 54 6209 20806 55 13915 4079 56 2441013196 57 13505 6117 58 9869 8220 59 1570 6044 60 25780 17387 61 2067124913 62 24558 20591 63 12402 3702 64 8314 1357 65 20071 14616 66 170143688 67 19837 946 68 15195 12136 69 7758 22808 70 3564 2925 71 3434 7769

[0084] TABLE 9 Address of Parity Bit Accumulators (Rate 8/9) 0 6235 28483222 1 5800 3492 5348 2 2757 927 90 3 6961 4516 4739 4 1172 3237 6264 51927 2425 3683 6 3714 6309 2495 7 3070 6342 7154 8 2428 613 3761 9 2906264 5927 10 1716 1950 4273 11 4613 6179 3491 12 4865 3286 6005 13 13435923 3529 14 4589 4035 2132 15 1579 3920 6737 16 1644 1191 5998 17 14822381 4620 18 6791 6014 6596 19 2738 5918 3786 0 5156 6166 1 1504 4356 2130 1904 3 6027 3187 4 6718 759 5 6240 2870 6 2343 1311 7 1039 5465 86617 2513 9 1588 5222 10 6561 535 11 4765 2054 12 5966 6892 13 1969 386914 3571 2420 15 4632 981 16 3215 4163 17 973 3117 18 3802 6198 19 37943948 0 3196 6126 1 573 1909 2 850 4034 3 5622 1601 4 6005 524 5 52515783 6 172 2032 7 1875 2475 8 497 1291 9 2566 3430 10 1249 740 11 29441948 12 6528 2899 13 2243 3616 14 867 3733 15 1374 4702 16 4698 2285 174760 3917 18 1859 4058 19 6141 3527 0 2148 5066 1 1306 145 2 2319 871 33463 1061 4 5554 6647 5 5837 339 6 5821 4932 7 6356 4756 8 3930 418 9211 3094 10 1007 4928 11 3584 1235 12 6982 2869 13 1612 1013 14 953 496415 4555 4410 16 4925 4842 17 5778 600 18 6509 2417 19 1260 4903 0 33693031 1 3557 3224 2 3028 583 3 3258 440 4 6226 6655 5 4895 1094 6 14816847 7 4433 1932 8 2107 1649 9 2119 2065 10 4003 6388 11 6720 3622 123694 4521 13 1164 7050 14 1965 3613 15 4331 66 16 2970 1796 17 4652 321818 1762 4777 19 5736 1399 0 970 2572 1 2062 6599 2 4597 4870 3 1228 69134 4159 1037 5 2916 2362 6 395 1226 7 6911 4548 8 4618 2241 9 4120 428010 5825 474 11 2154 5558 12 3793 5471 13 5707 1595 14 1403 325 15 66015183 16 6369 4569 17 4846 896 18 7092 6184 19 6764 7127 0 6358 1951 13117 6960 2 2710 7062 3 1133 3604 4 3694 657 5 1355 110 6 3329 6736 72505 3407 8 2462 4806 9 4216 214 10 5348 5619 11 6627 6243 12 2644 507313 4212 5088 14 3463 3889 15 5306 478 16 4320 6121 17 3961 1125 18 56991195 19 6511 792 0 3934 2778 1 3238 6587 2 1111 6596 3 1457 6226 4 14463885 5 3907 4043 6 6839 2873 7 1733 5615 8 5202 4269 9 3024 4722 10 54456372 11 370 1828 12 4695 1600 13 680 2074 14 1801 6690 15 2669 1377 162463 1681 17 5972 5171 18 5728 4284 19 1696 1459

[0085] TABLE 10 Address of Parity Bit Accumulators (Rate 9/10) 0 56112563 2900 1 5220 3143 4813 2 2481 834 81 3 6265 4064 4265 4 1055 29145638 5 1734 2182 3315 6 3342 5678 2246 7 2185 552 3385 8 2615 236 5334 91546 1755 3846 10 4154 5561 3142 11 4382 2957 5400 12 1209 5329 3179 131421 3528 6063 14 1480 1072 5398 15 3843 1777 4369 16 1334 2145 4163 172368 5055 260 0 6118 5405 1 2994 4370 2 3405 1669 3 4640 5550 4 13543921 5 117 1713 6 5425 2866 7 6047 683 8 5616 2582 9 2108 1179 10 9334921 11 5953 2261 12 1430 4699 13 5905 480 14 4289 1846 15 5374 6208 161775 3476 17 3216 2178 0 4165 884 1 2896 3744 2 874 2801 3 3423 5579 43404 3552 5 2876 5515 6 516 1719 7 765 3631 8 5059 1441 9 5629 598 105405 473 11 4724 5210 12 155 1832 13 1689 2229 14 449 1164 15 2308 308816 1122 669 17 2268 5758 0 5878 2609 1 782 3359 2 1231 4231 3 4225 20524 4286 3517 5 5531 3184 6 1935 4560 7 1174 131 8 3115 956 9 3129 1088 105238 4440 11 5722 4280 12 3540 375 13 191 2782 14 906 4432 15 3225 111116 6296 2583 17 1457 903 0 855 4475 1 4097 3970 2 4433 4361 3 5198 541 41146 4426 5 3202 2902 6 2724 525 7 1083 4124 8 2326 6003 9 5605 5990 104376 1579 11 4407 984 12 1332 6163 13 5359 3975 14 1907 1854 15 36015748 16 6056 3266 17 3322 4085 0 1768 3244 1 2149 144 2 1589 4291 3 51541252 4 1855 5939 5 4820 2706 6 1475 3360 7 4266 693 8 4156 2018 9 2103752 10 3710 3853 11 5123 931 12 6146 3323 13 1939 5002 14 5140 1437 151263 293 16 5949 4665 17 4548 6380 0 3171 4690 1 5204 2114 2 6384 5565 35722 1757 4 2805 6264 5 1202 2616 6 1018 3244 7 4018 5289 8 2257 3067 92483 3073 10 1196 5329 11 649 3918 12 3791 4581 13 5028 3803 14 31193506 15 4779 431 16 3888 5510 17 4387 4084 0 5836 1692 1 5126 1078 25721 6165 3 3540 2499 4 2225 6348 5 1044 1484 6 6323 4042 7 1313 5603 81303 3496 9 3516 3639 10 5161 2293 11 4682 3845 12 3045 643 13 2818 261614 3267 649 15 6236 593 16 646 2948 17 4213 1442 0 5779 1596 1 2403 12372 2217 1514 3 5609 716 4 5155 3858 5 1517 1312 6 2554 3158 7 5280 2643 84990 1353 9 5648 1170 10 1152 4366 11 3561 5368 12 3581 1411 13 56474661 14 1542 5401 15 5078 2687 16 316 1755 17 3392 1991

[0086] As regards the BCH encoder 211, the BCH code parameters areenumerated in Table 11. TABLE 11 BCH Uncoded Block BCH Error LDPC CodeLength BCH Coded Block Correction Rate k_(bch) Length n_(bch) (bits) 1/232208 32400 12 2/3 43040 43200 10 3/4 48408 48600 12 4/5 51648 51840 125/6 53840 54000 10 3/5 38688 38880 12 8/9 57472 57600 8  9/10 5819258320 8

[0087] It is noted that in the above table, n_(bch)=k_(ldpc).

[0088] The generator polynomial of the t error correcting BCH encoder211 is obtained by multiplying the first t polynomials in the followinglist of Table 12: TABLE 12 g1(x) 1 + x² + x³ + x⁵ + x¹⁶ g2(x) 1 + x +x⁴ + x⁵ + x⁶ + x⁸ + x¹⁶ g3(x) 1 + x² + x³ + x⁴ + x⁵ + x⁷ + x⁸ + x⁹ +x¹⁰ + x¹¹ + x¹⁶ g4(x) 1 + x² + x⁴ + x⁶ + x⁹ + x¹¹ + x¹² + x¹⁴ + x¹⁶g5(x) 1 + x + x² + x³ + x⁵ + x⁸ + x⁹ + x¹⁰ + x¹¹ + x¹² + x¹⁶ g6(x) 1 +x² + x⁴ + x⁵ + x⁷ + x⁸ + x⁹ + x¹⁰ + x¹² + x¹³ + x¹⁴ + x¹⁵ + x¹⁶ g7(x)1 + x² + x⁵ + x⁶ + x⁸ + x⁹ + x¹⁰ + x¹¹ + x¹³ + x¹⁵ + x¹⁶ g8(x) 1 + x +x² + x⁵ + x⁶ + x⁸ + x⁹ + x¹² + x¹³ + x¹⁴ + x¹⁶ g9(x) 1 + x⁵ + x⁷ + x⁹ +x¹⁰ + x¹¹ + x¹⁶ g10(x) 1 + x + x² + x⁵ + x⁷ + x⁸ + x¹⁰ + x¹² + x¹³ +x¹⁴ + x¹⁶ g11(x) 1 + x² + x³ + x⁵ + x⁹ + x¹¹ + x¹² + x¹³ + x¹⁶ g12(x)1 + x + x⁵ + x⁶ + x⁷ + x⁹ + x¹¹ + x¹² + x¹⁶

[0089] BCH encoding of information bits m=(m_(k) _(bch) ⁻¹, m_(k) _(bch)⁻², . . . , m₁, m₀) onto a codeword c=(m_(k) _(bch) ⁻¹, m_(k) _(bch) ⁻²,. . . , m₁, m₀, d_(n) _(bch) _(−k) _(bch) ⁻¹, d_(n) _(bch) _(−k) _(bch)⁻², . . . , d₁, d₀) is achieved as follows. The message polynomialm(x)=m_(k) _(bch) ⁻¹x^(k) ^(_(bch)) ⁻¹+m_(k) _(bch) ⁻²x^(k) ^(_(bch))⁻²+ . . . +m₁x+m₀ is multiplied by x^(n) ^(_(bch)) ^(−k) ^(_(bch)) .Next, x^(n) ^(_(bch)) ^(−k) ^(_(bch)) m(x) divided by g(x). Withd(x)=d_(n) _(bch) _(−k) _(bch) ⁻¹x^(n) ^(_(bch)) ^(−k) ^(_(bch)) ⁻¹+ . .. +d₁x+d₀ as the remainder, the codeword polynomial is set as follows:c(x)=x^(n) ^(_(bch)) ^(−k) ^(_(bch)) m(x)+d(x).

[0090] The above LDPC codes, in an exemplary embodiment, can be used tovariety of digital video applications, such as MPEG (Motion PicturesExpert Group) packet transmission.

[0091]FIG. 3 is a diagram of an exemplary receiver in the system ofFIG. 1. At the receiving side, a receiver 300 includes a demodulator 301that performs demodulation of received signals from transmitter 200.These signals are received at a receive antenna 303 for demodulation.After demodulation, the received signals are forwarded to a decoder 305,which attempts to reconstruct the original source messages by generatingmessages, X′, in conjunction with a bit metric generator 307. Withnon-Gray mapping, the bit metric generator 307 exchanges probabilityinformation with the decoder 305 back and forth (iteratively) during thedecoding process, which is detailed in FIG. 10. Alternatively, if Graymapping is used (according to one embodiment of the present invention),one pass of the bit metric generator is sufficient, in which furtherattempts of bit metric generation after each LDPC decoder iteration arelikely to yield limited performance improvement; this approach is morefully described with respect to FIG. 11. To appreciate the advantagesoffered by the present invention, it is instructive to examine how LDPCcodes are generated, as discussed in FIG. 4.

[0092]FIG. 4 is a diagram of a sparse parity check matrix, in accordancewith an embodiment of the present invention. LDPC codes are long, linearblock codes with sparse parity check matrix H_((n−k)xn). Typically theblock length, n, ranges from thousands to tens of thousands of bits. Forexample, a parity check matrix for an LDPC code of length n=8 and rate ½is shown in FIG. 4. The same code can be equivalently represented by thebipartite graph, per FIG. 5.

[0093]FIG. 5 is a diagram of a bipartite graph of an LDPC code of thematrix of FIG. 4. Parity check equations imply that for each check node,the sum (over GF (Galois Field)(2)) of all adjacent bit nodes is equalto zero. As seen in the figure, bit nodes occupy the left side of thegraph and are associated with one or more check nodes, according to apredetermined relationship. For example, corresponding to check node m₁,the following expression exists n₁+n₄+n₅+n₈+0 with respect to the bitnodes.

[0094] Returning the receiver 303, the LDPC decoder 305 is considered amessage passing decoder, whereby the decoder 305 aims to find the valuesof bit nodes. To accomplish this task, bit nodes and check nodesiteratively communicate with each other. The nature of thiscommunication is described below.

[0095] From check nodes to bit nodes, each check node provides to anadjacent bit node an estimate (“opinion”) regarding the value of thatbit node based on the information coming from other adjacent bit nodes.For instance, in the above example if the sum of n₄, n₅ and n₈ “lookslike” 0 to m₁, then m₁ would indicate to n₁ that the value of n₁ isbelieved to be 0 (since n₁+n₄+n₅+n₈=0); otherwise m₁ indicate to n₁ thatthe value of n₁ is believed to be 1. Additionally, for soft decisiondecoding, a reliability measure is added.

[0096] From bit nodes to check nodes, each bit node relays to anadjacent check node an estimate about its own value based on thefeedback coming from its other adjacent check nodes. In the aboveexample n, has only two adjacent check nodes m₁ and m₃. If the feedbackcoming from m₃ to n₁ indicates that the value of n, is probably 0, thenn₁ would notify m₁ that an estimate of n₁'s own value is 0. For the casein which the bit node has more than two adjacent check nodes, the bitnode performs a majority vote (soft decision) on the feedback comingfrom its other adjacent check nodes before reporting that decision tothe check node it communicates. The above process is repeated until allbit nodes are considered to be correct (i.e., all parity check equationsare satisfied) or until a predetermined maximum number of iterations isreached, whereby a decoding failure is declared.

[0097]FIG. 6 is a diagram of a sub-matrix of a sparse parity checkmatrix, wherein the sub-matrix contains parity check values restrictedto the lower triangular region, according to an embodiment of thepresent invention. As described previously, the encoder 203 (of FIG. 2)can employ a simple encoding technique by restricting the values of thelower triangular area of the parity check matrix. According to anembodiment of the present invention, the restriction imposed on theparity check matrix is of the form:

H _((n−k)xn) =[A _((n−k)xk) B _((n−k)x(n−k))]

[0098] , where B is lower triangular.

[0099] Any information block i=(i₀, i₁, . . . , i_(k−1)) is encoded to acodeword c=(i₀, i₁, . . . , i_(k−1), p₀, p₁, . . . p_(n−k−1)) usingHc^(T)=0, and recursively solving for parity bits; for example,

a ₀₀ i ₀ +a ₀₁ i ₁ + . . . +a _(0,k−1) i _(k−1) +p ₀=0

Solve p ₀

a ₁₀ i ₀ +a ₁₁ i ₁ + . . . +a _(1,k−1) i _(k−1) +b ₁₀ p ₀ +p ₁=0

[0100] and similarly for p₂, p₃, . . . , p_(n−k−1).

[0101]FIG. 7 is a graph showing performance between codes utilizingunrestricted parity check matrix (H matrix) versus restricted H matrixof FIG. 6. The graph shows the performance comparison between two LDPCcodes: one with a general parity check matrix and the other with aparity check matrix restricted to be lower triangular to simplifyencoding. The modulation scheme, for this simulation, is 8-PSK. Theperformance loss is within 0.1 dB. Therefore, the performance loss isnegligible based on the restriction of the lower triangular H matrices,while the gain in simplicity of the encoding technique is significant.Accordingly, any parity check matrix that is equivalent to a lowertriangular or upper triangular under row and/or column permutation canbe utilized for the same purpose.

[0102]FIGS. 8A and 8B are, respectively, a diagram of a non-Gray 8-PSKmodulation scheme, and a Gray 8-PSK modulation, each of which can beused in the system of FIG. 1. The non-Gray 8-PSK scheme of FIG. 8A canbe utilized in the receiver of FIG. 3 to provide a system that requiresvery low Frame Erasure Rate (FER). This requirement can also besatisfied by using a Gray 8-PSK scheme, as shown in FIG. 8B, inconjunction with an outer code, such as Bose, Chaudhuri, and Hocquenghem(BCH), Hamming, or Reed-Solomon (RS) code.

[0103] Under this scheme, there is no need to iterate between the LDPCdecoder 305 (FIG. 3) and the bit metric generator 307, which may employ8-PSK modulation. In the absence of an outer code, the LDPC decoder 305using Gray labeling exhibit an earlier error floor, as shown in FIG. 9below.

[0104]FIG. 8C shows a diagram of a process for bit labeling for a higherorder signal constellation, according to an embodiment of the presentinvention. A codeword is output from the LDPC encoder 203 (FIGS. 2A and2B), and is mapped to a constellation point in a higher order signalconstellation (as shown in FIGS. 8D and 8F), per steps 801, 803. Thismapping is not performed sequentially as in traditional systems, butinstead executed on a non-sequential basis, such as interleaving. Such amapping is further detailed below with respect to FIG. 8F. The modulator205 then modulates, as in step 805, the signal based on the mapping. Themodulated signal is thereafter transmitted (step 807).

[0105]FIG. 8D shows a diagram of exemplary 16-APSK (Amplitude PhaseShift Keying) constellations. Constellations A and B are 16-APSKconstellations. The only difference between the two constellations A andB is that the inner circle symbols of Constellation A are rotated 15degrees counterclockwise with respect to the inner circle symbols ofConstellation B, such that inner circle symbols fall between the outercircle symbols to maximize inter-symbol distances. Therefore,intuitively Constellation A is more attractive if the Forward ErrorCorrection (FEC) decoder 305 used a symbolwise decoding algorithm. Onthe other hand, given the multiplicity of code rates and differentconstellations, using an FEC code tailored towards bitwise decoding ismore flexible. In such a case, it is not apparent which constellationswould perform better, in that while Constellation A maximizes symbolwisedistances, Constellation B is more “Gray-coding friendly.” AWGN(Additive White Gaussian Noise) simulations, with code rate ¾, wereperformed (the results of which are shown in FIG. 8E) that with bitwisedecoding, Constellation B performs slightly better.

[0106]FIG. 8F is a diagram of constellations for Quadrature Phase ShiftKeying (QPSK), 8-PSK, 16-APSK and 32-APSK symbols, in accordance with anembodiment of the present invention;

[0107]FIG. 8F show symmetric constellations for QPSK, 8-PSK, 16-APSK and32-APSK symbols, respectively. With QSPK, two LDPC coded bits from theLDPC encoder 203 are mapped to a QPSK symbol. That is, bits 2 i and 2i+1 determines the i^(th) QPSK symbol, where i=0, 1, 2 . . . , N/2−1,and N is the coded LDPC block size. For 8-PSK, bits N/3+i, 2N/3+i and idetermine the i^(th) 8-PSK symbol, where i=0, 1, 2, . . . , N/3−1. For16-APSK, bits N/2+2i, 2i, N/2+2i+1 and 2 i+1 specify the i^(th) 16-APSKsymbol, where i=0, 1, 2, . . . , N/4−1. Further, for 32-APSK, bitsN/5+i, 2N/5+i, 4N/5+i, 3N/5+i and i determine the i^(th) symbol, wherei=0, 1, 2, . . . , N/5−1.

[0108] Alternatively, 8-PSK, 16-APSK and 32-APSK constellation labelingcan be chosen as shown in FIG. 8G. With this labeling, N LDPC encodedbits are first passed through a bit interleaver. The bit interleavingtable, in an exemplary embodiment, is a two-dimensional array with N/3rows and 3 columns for 8-PSK, N/4 rows and 4 columns for 16-APSK and N/5rows and 5 columns for 32-APSK. The LDPC encoded bits are written to theinterleaver table column by column, and read out row by row. It is notedthat for the case of 8-PSK and 32-APSK, this row/column bit interleaverstrategy with labeling as shown in FIG. 8G, is exactly equivalent to thebit interleaving strategy described above with respect to the labelingshown in FIG. 8F. For the case of 16-APSK, these two strategies arefunctionally equivalent; that is, they exhibit the same performance onan AWGN channel.

[0109]FIG. 8H illustrates the simulation results (on AWGN Channel) ofthe above symbol constellations. Table 13 summarizes expectedperformance at PER=10⁻⁶ and distance from constrained capacity. TABLE 13Distance to Code Es/No (dB) Capacity (dB) 2/3, 8-PSK 6.59 0.873 3/4,8-PSK 7.88 0.690 5/6, 8-PSK 9.34 0.659 8/9, 8-PSK 10.65 0.750 9/10,8-PSK 10.95 0.750 1/2, QPSK 0.99 0.846 3/5, QPSK 2.20 0.750 2/3, QPSK3.07 0.760 3/4, QPSK 4.02 0.677 4/5, QPSK 4.66 0.627 5/6, QPSK 5.150.600 7/8, QPSK 5.93 0.698 8/9, QPSK 6.17 0.681 9/10, QPSK 6.39 0.6873/4, 16-APSK 10.19 0.890 4/5, 16-APSK 11.0 0.850 5/6, 16-APSK 11.580.800 7/8, 16-APSK 12.54 0.890 4/5, 32-APSK 13.63 1.100 5/6, 32-APSK14.25 1.050 8/9, 32-APSK 15.65 1.150

[0110]FIG. 9 is a graph showing performance between codes utilizing Graylabeling versus non-Gray labeling of FIGS. 8A and 8B. The error floorstems from the fact that assuming correct feedback from LDPC decoder305, regeneration of 8-PSK bit metrics is more accurate with non-Graylabeling since the two 8-PSK symbols with known two bits are furtherapart with non-Gray labeling. This can be equivalently seen as operatingat higher Signal-to-Noise Ratio (SNR). Therefore, even though errorasymptotes of the same LDPC code using Gray or non-Gray labeling havethe same slope (i.e., parallel to each other), the one with non-Graylabeling passes through lower FER at any SNR.

[0111] On the other hand, for systems that do not require very low FER,Gray labeling without any iteration between LDPC decoder 305 and 8-PSKbit metric generator 307 may be more suitable because re-generating8-PSK bit metrics before every LDPC decoder iteration causes additionalcomplexity. Moreover, when Gray labeling is used, re-generating 8-PSKbit metrics before every LDPC decoder iteration yields only very slightperformance improvement. As mentioned previously, Gray labeling withoutiteration may be used for systems that require very low FER, provided anouter code is implemented.

[0112] The choice between Gray labeling and non-Gray labeling dependsalso on the characteristics of the LDPC code. Typically, the higher bitor check node degrees, the better it is for Gray labeling, because forhigher node degrees, the initial feedback from LDPC decoder 305 to 8-PSK(or similar higher order modulation) bit metric generator 307deteriorates more with non-Gray labeling.

[0113] When 8-PSK (or similar higher order) modulation is utilized witha binary decoder, it is recognized that the three (or more) bits of asymbol are not received “equally noisy”. For example with Gray 8-PSKlabeling, the third bit of a symbol is considered more noisy to thedecoder than the other two bits. Therefore, the LDPC code design doesnot assign a small number of edges to those bit nodes represented by“more noisy” third bits of 8-PSK symbol so that those bits are notpenalized twice.

[0114]FIG. 10 is a flow chart of the operation of the LDPC decoder usingnon-Gray mapping, according to an embodiment of the present invention.Under this approach, the LDPC decoder and bit metric generator iterateone after the other. In this example, 8-PSK modulation is utilized;however, the same principles apply to other higher modulation schemes aswell. Under this scenario, it is assumed that the demodulator 301outputs a distance vector, d, denoting the distances between receivednoisy symbol points and 8-PSK symbol points to the bit metric generator307, whereby the vector components are as follows:$d_{i} = {- \frac{E_{s}}{N_{0}}}$

 {(r _(x) −s _(i,x))²+(r _(y) −s _(i,y))²} i=0, 1, . . . 7.

[0115] The 8-PSK bit metric generator 307 communicates with the LDPCdecoder 305 to exchange a priori probability information and aposteriori probability information, which respectively are representedas u, and a. That is, the vectors u and a respectively represent apriori and a posteriori probabilities of log likelihood ratios of codedbits.

[0116] The 8-PSK bit metric generator 307 generates the a priorilikelihood ratios for each group of three bits as follows. First,extrinsic information on coded bits is obtained:

e _(j) =a _(j) −u _(j) j=0, 1, 2.

[0117] Next, 8-PSK symbol probabilities, Pi i=0, 1, . . . , 7, aredetermined.

[0118] *y_(j)=−ƒ(0, e_(j)) j=0,1,2 where ƒ(a, b)=max(a,b)+LUT_(ƒ)(a, b)with LUT_(ƒ)(a,b)=ln(l+e^(−|a−b|))

[0119] *x_(j)=y_(j)+e_(j) j=0, 1, 2

[0120] *p₀=x₀+x₁+x₂ p₄=y₀+x₁+x₂

[0121] p₁=x₀+x₁+y₂ p₅=y₀+x₁+y₂

[0122] p₂=x₀+y₁+x₂ p₆=y₀+y₁+x₂

[0123] p₃=x₀+y₁+y₂ p₇=y₀+y₁+y₂

[0124] Next, the bit metric generator 307 determines a priori loglikelihood ratios of the coded bits as input to LDPC decoder 305, asfollows:

u ₀=ƒ(d ₀ +p ₀ ,d ₁ +p ₁ ,d ₂ +p ₂ ,d ₃ +p ₃)−ƒ(d ₄ +p ₄ ,d ₅ +p ₅ ,d ₆+p ₆ ,d ₇ +p ₇)−e ₀

u ₁=ƒ(d ₀ +p ₀ ,d ₁ +p ₁ ,d ₄ +p ₄ ,d ₅ +p ₅)−ƒ(d ₂ +p ₂ ,d ₃ +p ₃ ,d₆+p₆ ,d ₇ +p ₇)−e ₁

u ₂=ƒ(d ₀ +p ₀ ,d ₂ +p ₂ ,d ₄ +p ₄ ,d ₆ +p ₆)−ƒ(d ₁ +p ₁ ,d ₃ +p ₃ ,d ₅+p ₅ ,d ₇ +p ₇)−e ₂

[0125] It is noted that the function ƒ(.) with more than two variablescan be evaluated recursively; e.g. ƒ(a,b,c)=ƒ(ƒ(a,b),c).

[0126] The operation of the LDPC decoder 305 utilizing non-Gray mappingis now described. In step 1001, the LDPC decoder 305 initializes loglikelihood ratios of coded bits, v, before the first iteration accordingto the following (and as shown in FIG. 12A):

v_(n→k) _(i) =u_(n), n=0, 1, . . . , N−1, i=1, 2, . . . , deg(bit noden)

[0127] Here, v_(n→k) _(i) denotes the message that goes from bit node nto its adjacent check node k_(i), u_(n), denotes the demodulator outputfor the bit n and N is the codeword size.

[0128] In step 1003, a check node, k, is updated, whereby the input vyields the output w. As seen in FIG. 12B, the incoming messages to thecheck node k from its d_(c) adjacent bit nodes are denoted by v_(n) ₁_(→k), v_(n) _(2→k) , . . . , v_(n) _(dc) _(→k). The goal is to computethe outgoing messages from the check node k back to d_(c) adjacent bitnodes. These messages are denoted by w_(k→n) ₁ , w_(k→n) ₂ , . . . ,w_(k→n) _(dc) , where

w _(k→n) ₁ =g(v _(n) ₁ _(→k) ,v _(n) ₂ _(→k) , . . . ,v _(n) _(i−1)_(→k) , . . . ,v _(n) _(dc) _(→k)).

[0129] The function g( ) is defined as follows:

g(a,b)=sign(a)×sign(b)×{min(|a|,|b|)}+LUT _(g)(a,b),

[0130] where LUT_(g)(a, b)=ln(l+e^(−|a+b|))−ln(l+e^(−|a−b|)) Similar tofunction ƒ, function g with more than two variables can be evaluatedrecursively.

[0131] Next, the decoder 305, per step 1205, outputs a posterioriprobability information (FIG. 12C), such that:$a_{n} = {u_{n} + {\sum\limits_{j}{w_{k_{j}\rightarrow n}.}}}$

[0132] Per step 1007, it is determined whether all the parity checkequations are satisfied. If these parity check equations are notsatisfied, then the decoder 305, as in step 1009, re-derives 8-PSK bitmetrics and channel input u_(n). Next, the bit node is updated, as instep 1011. As shown in FIG. 14C, the incoming messages to the bit node nfrom its d, adjacent check nodes are denoted by w_(k) ₁ _(→n), w_(k) ₂_(→n), . . . , w_(k) _(dv) _(→n) The outgoing messages from the bit noden are computed back to d_(v) adjacent check nodes; such messages aredenoted by v_(n→k) ₁ , v_(n→k) ₂ , . . . , v_(n→k) _(dv) , and computedas follows:$v_{n\rightarrow k_{i}} = {u_{n} + {\sum\limits_{j \neq i}w_{k_{j}\rightarrow n}}}$

[0133] In step 1013, the decoder 305 outputs the hard decision (in thecase that all parity check equations are satisfied):${\hat{c}}_{n} = \left\{ \begin{matrix}{0,{a_{n} \geq 0}} \\{1,{a_{n} < 0}}\end{matrix} \right.$

[0134] Stop if Hĉ^(T)=0

[0135] The above approach is appropriate when non-Gray labeling isutilized. However, when Gray labeling is implemented, the process ofFIG. 11 is executed.

[0136]FIG. 11 is a flow chart of the operation of the LDPC decoder ofFIG. 3 using Gray mapping, according to an embodiment of the presentinvention. When Gray labeling is used, bit metrics are advantageouslygenerated only once before the LDPC decoder, as re-generating bitmetrics after every LDPC decoder iteration may yield nominal performanceimprovement. As with steps 1001 and 1003 of FIG. 10, initialization ofthe log likelihood ratios of coded bits, v, are performed, and the checknode is updated, per steps 1101 and 1103. Next, the bit node n isupdated, as in step 1105. Thereafter, the decoder outputs the aposteriori probability information (step 1107). In step 1109, adetermination is made whether all of the parity check equations aresatisfied; if so, the decoder outputs the hard decision (step 1111).Otherwise, steps 1103-1107 are repeated.

[0137]FIG. 13A is a flowchart of process for computing outgoing messagesbetween the check nodes and the bit nodes using a forward-backwardapproach, according to an embodiment of the present invention. For acheck node with d_(c) adjacent edges, the computation of d_(c)(d_(c)−1)and numerous g(.,.) functions are performed. However, theforward-backward approach reduces the complexity of the computation to3(d_(c) −2), in which d _(c)−1 variables are stored.

[0138] Referring to FIG. 12B, the incoming messages to the check node kfrom d, adjacent bit nodes are denoted by v_(n) _(1→k, v) _(n) ₂ _(→k),. . . , v_(n) _(de) _(→k). It is desired that the outgoing messages arecomputed from the check node k back to d_(c) adjacent bit nodes; theseoutgoing messages are denoted by w_(k→n) ₁ , w_(k→n) ₂ , . . . , w_(k→n)_(dc) .

[0139] Under the forward-backward approach to computing these outgoingmessages, forward variables, ƒ₁, ƒ₂, . . . , ƒ_(dc), are defined asfollows:

ƒ₁=v_(1→k)

ƒ₂ =g(ƒ₁ ,v _(2→k))

ƒ₃ =g(ƒ₂ ,v _(3→k))

: : :

ƒ_(dc) =g(ƒ_(dc−1) ,v _(dc→k))

[0140] In step 1301, these forward variables are computed, and stored,per step 1303.

[0141] Similarly, backward variables, b₁, b₂, . . . , b_(dc), aredefined by the following:

b _(dc) =v _(dc→k)

b _(dc−1) =g(b _(dc) ,v _(dc−1→k))

: : :

b ₁ =g(b ₂ ,v _(1→k))

[0142] In step 1305, these backward variables are then computed.Thereafter, the outgoing messages are computed, as in step 1307, basedon the stored forward variables and the computed backward variables. Theoutgoing messages are computed as follows:

w_(k→1)=b₂

w _(k→i) =g(f _(i−1) ,b _(i+1)) i=2, 3, . . . , d_(c)−1

w _(k→dc)=f_(dc−1)

[0143] Under this approach, only the forward variables, f₂, f₃, . . . ,f_(dc), are required to be stored. As the backward variables b_(i) arecomputed, the outgoing messages, w_(k→i), are simultaneously computed,thereby negating the need for storage of the backward variables.

[0144] The computation load can be further enhance by a parallelapproach, as next discussed.

[0145]FIG. 13B is a flowchart of process for computing outgoing messagesbetween the check nodes and the bit nodes using a parallel approach,according to an embodiment of the present invention. For a check node kwith inputs v_(n) _(1→k) , v_(n) _(2→k) , . . . , v_(n) _(dc) _(→k) fromd_(c) adjacent bit nodes, the following parameter is computed, as instep 1311:

γ_(k) =g(v _(n) ₁ _(→k) ,v _(n) ₂ _(→k) , . . . ,v _(n) _(dc) _(→k)).

[0146] It is noted that the g(.,.) function can also be expressed asfollows:${g\left( {a,b} \right)} = {\ln {\frac{1 + ^{a + b}}{^{a} + ^{b}}.}}$

[0147] Exploiting the recursive nature of the g(.,.) function, thefollowing expression results:$\gamma_{k} = {{\ln \frac{1 + ^{{g{({{v_{{n_{1}\rightarrow k},\quad \ldots \quad,}v_{n_{i - 1}\rightarrow k}},v_{n_{i + 1}\rightarrow k},\quad \ldots \quad,v_{n_{d\quad c}\rightarrow k}})}} + v_{n_{i}\rightarrow k}}}{^{{g{({{v_{{n_{j}\rightarrow k},\quad \ldots \quad,}v_{n_{i - j}\rightarrow k}},v_{n_{i + 1}\rightarrow k},\quad \ldots \quad,v_{n_{d\quad c}\rightarrow k}})}} + ^{v_{n_{i}\rightarrow k}}}}} = {\ln \frac{1 + ^{w_{k\rightarrow n_{i}} + v_{n_{i}\rightarrow k}}}{^{w_{k\rightarrow n_{i}}} + ^{v_{n_{i}\rightarrow k}}}}}$

[0148] Accordingly, w_(k→n) _(i) can be solved in the following manner:$w_{k\rightarrow n_{i}} = {{\ln \frac{^{v_{n_{i}\rightarrow k} + \gamma_{k}} - 1}{^{v_{n_{i}\rightarrow k} - \gamma_{k}} - 1}} - \gamma_{k}}$

[0149] The ln(.) term of the above equation can be obtained using alook-up table LUT_(x) that represents the function ln |e^(x)−1| (step1313). Unlike the other look-up tables LUT_(f) or LUT_(g), the tableLUT_(x) would likely requires as many entries as the number ofquantization levels. Once γ_(k) is obtained, the calculation of w_(k→n)_(i) for all n_(i) can occur in parallel using the above equation, perstep 1315.

[0150] The computational latency of γ_(k) is advantageously 1og₂(d_(c)).

[0151] FIGS. 14A-14C are graphs showing simulation results of LDPC codesgenerated in accordance with various embodiments of the presentinvention. In particular, FIGS. 14A-14C show the performance of LDPCcodes with higher order modulation and code rates of ¾ (QPSK, 1.485bits/symbol), ⅔ (8-PSK, 1.980 bits/symbol), and ⅚ (8-PSK, 2.474bits/symbol).

[0152] Two general approaches exist to realize the interconnectionsbetween check nodes and bit nodes: (1) a fully parallel approach, and(2) a partially parallel approach. In fully parallel architecture, allof the nodes and their interconnections are physically implemented. Theadvantage of this architecture is speed.

[0153] The fully parallel architecture, however, may involve greatercomplexity in realizing all of the nodes and their connections.Therefore with fully parallel architecture, a smaller block size may berequired to reduce the complexity. In that case, for the same clockfrequency, a proportional reduction in throughput and some degradationin FER versus Es/No performance may result.

[0154] The second approach to implementing LDPC codes is to physicallyrealize only a subset of the total number of the nodes and use onlythese limited number of “physical” nodes to process all of the“functional” nodes of the code. Even though the LDPC decoder operationscan be made extremely simple and can be performed in parallel, thefurther challenge in the design is how the communication is establishedbetween “randomly” distributed bit nodes and check nodes. The decoder305 (of FIG. 3), according to one embodiment of the present invention,addresses this problem by accessing memory in a structured way, as torealize a seemingly random code. This approach is explained with respectto, FIGS. 15A and 15B.

[0155]FIGS. 15A and 15B are diagrams of the top edge and bottom edge,respectively, of memory organized to support structured access as torealize randomness in LDPC coding, according to an embodiment of thepresent invention. Structured access can be achieved withoutcompromising the performance of a truly random code by focusing on thegeneration of the parity check matrix. In general, a parity check matrixcan be specified by the connections of the check nodes with the bitnodes. For example, the bit nodes can be divided into groups of a fixedsize, which for illustrative purposes is 392. Additionally, assuming thecheck nodes connected to the first bit node of degree 3, for instance,are numbered as a, b and c, then the check nodes connected to the secondbit node are numbered as a+p, b+p and c+p, the check nodes connected tothe third bit node are numbered as a+2p, b+2p and c+2p etc.; wherep=(number of check nodes)/392. For the next group of 392 bit nodes, thecheck nodes connected to the first bit node are different from a, b, cso that with a suitable choice of p, all the check nodes have the samedegree. A random search is performed over the free constants such thatthe resulting LDPC code is cycle-4 and cycle-6 free. Because of thestructural characteristics of the parity check matrix of the presentinvention, the edge information can stored to permit concurrent accessto a group of relevant edge values during decoding.

[0156] In other words, the approach of the present invention facilitatesmemory access during check node and bit node processing. The values ofthe edges in the bipartite graph can be stored in a storage medium, suchas random access memory (RAM). It is noted that for a truly random LDPCcode during check node and bit node processing, the values of the edgeswould need to be accessed one by one in a random fashion. However, sucha conventional access scheme would be too slow for a high data rateapplication. The RAM of FIGS. 15A and 15B are organized in a manner,whereby a large group of relevant edges can be fetched in one clockcycle; accordingly, these values are placed “together” in memory,according to a predetermined scheme or arrangement. It is observed that,in actuality, even with a truly random code, for a group of check nodes(and respectively bit nodes), the relevant edges can be placed next toone another in RAM, but then the relevant edges adjacent to a group ofbit nodes (respectively check nodes) will be randomly scattered in RAM.Therefore, the “togetherness,” under the present invention, stems fromthe design of the parity check matrices themselves. That is, the checkmatrix design ensures that the relevant edges for a group of bit nodesand check nodes are simultaneously placed together in RAM.

[0157] As seen in FIGS. 15A and 15B, each box contains the value of anedge, which is multiple bits (e.g., 6). Edge RAM, according to oneembodiment of the present invention, is divided into two parts: top edgeRAM 1501 (FIG. 15A) and bottom edge RAM 1503 (FIG. 15B). Bottom edge RAMcontains the edges between bit nodes of degree 2, for example, and checknodes. Top edge RAM 1501 contains the edges between bit nodes of degreegreater than 2 and check nodes. Therefore, for every check node, 2adjacent edges are stored in the bottom RAM 1503, and the rest of theedges are stored in the top edge RAM 1501. For example, the size of thetop edge RAM 1501 and bottom edge RAM 1503 for various code rates aregiven in Table 14: TABLE 14 1/2 2/3 3/4 5/6 TopEdge 400 × 392 440 × 392504 × 392 520 × 392 RAM Bottom 160 × 392 110 × 392  72 × 392  52 × 392Edge RAM

[0158] Based on Table 14, an edge RAM of size 576×392 is sufficient tostore the edge metrics for all the code rates of ½, ⅔, ¾, and ⅚.

[0159] As noted, under this exemplary scenario, a group of 392 bit nodesand 392 check nodes are selected for processing at a time. For 392 checknode processing, q=d_(c)−2 consecutive rows are accessed from the topedge RAM 1501, and 2 consecutive rows from the bottom edge RAM 1503. Thevalue of dc depends on the specific code, for example d_(c)=7 for rate½, d_(c)=10 for rate ⅔, d_(c)=16 for rate ¾ and d_(c)=22 for rate ⅚ forthe above codes. Of course other values of dc for other codes arepossible. In this instance, q+2 is the degree of each check node.

[0160] For bit node processing, if the group of 392 bit nodes has degree2, their edges are located in 2 consecutive rows of the bottom edge RAM1503. If the bit nodes have degree d>2, their edges are located in somed rows of the top edge RAM 1501. The address of these d rows can bestored in non-volatile memory, such as Read-Only Memory (ROM). The edgesin one of the rows correspond to the first edges of 392 bit nodes, theedges in another row correspond to the second edges of 392 bit nodes,etc. Moreover for each row, the column index of the edge that belongs tothe first bit node in the group of 392 can also be stored in ROM. Theedges that correspond to the second, third, etc. bit nodes follow thestarting column index in a “wrapped around” fashion. For example, if thej^(th) edge in the row belongs to the first bit node, then the (j+1)stedge belongs to the second bit node, (j+2)nd edge belongs to the thirdbit node . . . , and (j−1)st edge belongs to the 392^(st) bit node.

[0161] With the organization shown in FIGS. 15A and 15B, speed of memoryaccess is greatly enhanced during LDPC coding.

[0162]FIG. 16 illustrates a computer system upon which an embodimentaccording to the present invention can be implemented. The computersystem 1600 includes a bus 1601 or other communication mechanism forcommunicating information, and a processor 1603 coupled to the bus 1601for processing information. The computer system 1600 also includes mainmemory 1605, such as a random access memory (RAM) or other dynamicstorage device, coupled to the bus 1601 for storing information andinstructions to be executed by the processor 1603. Main memory 1605 canalso be used for storing temporary variables or other intermediateinformation during execution of instructions to be executed by theprocessor 1603. The computer system 1600 further includes a read onlymemory (ROM) 1607 or other static storage device coupled to the bus 1601for storing static information and instructions for the processor 1603.A storage device 1609, such as a magnetic disk or optical disk, isadditionally coupled to the bus 1601 for storing information andinstructions.

[0163] The computer system 1600 may be coupled via the bus 1601 to adisplay 1611, such as a cathode ray tube (CRT), liquid crystal display,active matrix display, or plasma display, for displaying information toa computer user. An input device 1613, such as a keyboard includingalphanumeric and other keys, is coupled to the bus 1601 forcommunicating information and command selections to the processor 1603.Another type of user input device is cursor control 1615, such as amouse, a trackball, or cursor direction keys for communicating directioninformation and command selections to the processor 1603 and forcontrolling cursor movement on the display 1611.

[0164] According to one embodiment of the invention, generation of LDPCcodes is provided by the computer system 1600 in response to theprocessor 1603 executing an arrangement of instructions contained inmain memory 1605. Such instructions can be read into main memory 1605from another computer-readable medium, such as the storage device 1609.Execution of the arrangement of instructions contained in main memory1605 causes the processor 1603 to perform the process steps describedherein. One or more processors in a multi-processing arrangement mayalso be employed to execute the instructions contained in main memory1605. In alternative embodiments, hard-wired circuitry may be used inplace of or in combination with software instructions to implement theembodiment of the present invention. Thus, embodiments of the presentinvention are not limited to any specific combination of hardwarecircuitry and software.

[0165] The computer system 1600 also includes a communication interface1617 coupled to bus 1601. The communication interface 1617 provides atwo-way data communication coupling to a network link 1619 connected toa local network 1621. For example, the communication interface 1617 maybe a digital subscriber line (DSL) card or modem, an integrated servicesdigital network (ISDN) card, a cable modem, or a telephone modem toprovide a data communication connection to a corresponding type oftelephone line. As another example, communication interface 1617 may bea local area network (LAN) card (e.g. for Ethernet™ or an AsynchronousTransfer Model (ATM) network) to provide a data communication connectionto a compatible LAN. Wireless links can also be implemented. In any suchimplementation, communication interface 1617 sends and receiveselectrical, electromagnetic, or optical signals that carry digital datastreams representing various types of information. Further, thecommunication interface 1617 can include peripheral interface devices,such as a Universal Serial Bus (USB) interface, a PCMCIA (PersonalComputer Memory Card International Association) interface, etc.

[0166] The network link 1619 typically provides data communicationthrough one or more networks to other data devices. For example, thenetwork link 1619 may provide a connection through local network 1621 toa host computer 1623, which has connectivity to a network 1625 (e.g. awide area network (WAN) or the global packet data communication networknow commonly referred to as the “Internet”) or to data equipmentoperated by service provider. The local network 1621 and network 1625both use electrical, electromagnetic, or optical signals to conveyinformation and instructions. The signals through the various networksand the signals on network link 1619 and through communication interface1617, which communicate digital data with computer system 1600, areexemplary forms of carrier waves bearing the information andinstructions.

[0167] The computer system 1600 can send messages and receive data,including program code, through the network(s), network link 1619, andcommunication interface 1617. In the Internet example, a server (notshown) might transmit requested code belonging to an application programfor implementing an embodiment of the present invention through thenetwork 1625, local network 1621 and communication interface 1617. Theprocessor 1603 may execute the transmitted code while being receivedand/or store the code in storage device 169, or other non-volatilestorage for later execution. In this manner, computer system 1600 mayobtain application code in the form of a carrier wave.

[0168] The term “computer-readable medium” as used herein refers to anymedium that participates in providing instructions to the processor 1603for execution. Such a medium may take many forms, including but notlimited to non-volatile media, volatile media, and transmission media.Non-volatile media include, for example, optical or magnetic disks, suchas storage device 1609. Volatile media include dynamic memory, such asmain memory 1605. Transmission media include coaxial cables, copper wireand fiber optics, including the wires that comprise bus 1601.Transmission media can also take the form of acoustic, optical, orelectromagnetic waves, such as those generated during radio frequency(RF) and infrared (IR) data communications. Common forms ofcomputer-readable media include, for example, a floppy disk, a flexibledisk, hard disk, magnetic tape, any other magnetic medium, a CD-ROM,CDRW, DVD, any other optical medium, punch cards, paper tape, opticalmark sheets, any other physical medium with patterns of holes or otheroptically recognizable indicia, a RAM, a PROM, and EPROM, a FLASH-EPROM,any other memory chip or cartridge, a carrier wave, or any other mediumfrom which a computer can read.

[0169] Various forms of computer-readable media may be involved inproviding instructions to a processor for execution. For example, theinstructions for carrying out at least part of the present invention mayinitially be borne on a magnetic disk of a remote computer. In such ascenario, the remote computer loads the instructions into main memoryand sends the instructions over a telephone line using a modem. A modemof a local computer system receives the data on the telephone line anduses an infrared transmitter to convert the data to an infrared signaland transmit the infrared signal to a portable computing device, such asa personal digital assistance (PDA) and a laptop. An infrared detectoron the portable computing device receives the information andinstructions borne by the infrared signal and places the data on a bus.The bus conveys the data to main memory, from which a processorretrieves and executes the instructions. The instructions received bymain memory may optionally be stored on storage device either before orafter execution by processor.

[0170] Accordingly, the various embodiments of the present inventionprovide an approach is provided for bit labeling of a signalconstellation. An encoder, such as a Low Density Parity Check (LDPC)encoder, generates encoded signals by transforming an input message intoa codeword represented by a plurality of set of bits. These bits aremapped non-sequentially (e.g., interleaving) a higher orderconstellation (Quadrature Phase Shift Keying (QPSK), 8-PSK, 16-APSK(Amplitude Phase Shift Keying), 32-APSK, etc. The above arrangementadvantageously provides enhanced performance of the codes.

[0171] While the present invention has been described in connection witha number of embodiments and implementations, the present invention isnot so limited but covers various obvious modifications and equivalentarrangements, which fall within the purview of the appended claims.

What is claimed is:
 1. A method for transmitting encoded signals, themethod comprising: receiving one of a plurality of set of bits of acodeword from an encoder for transforming an input message into thecodeword; non-sequentially mapping the one set of bits into a higherorder constellation; and outputting a symbol of the higher orderconstellation corresponding to the one set of bits based on the mapping.2. A method according to claim 1, further comprising: writing N encodedbits to a block interleaver on a column by column basis; and reading outthe encoded bits on a row by row basis, wherein the block interleaverhas N/3 rows and 3 columns when the higher order modulation is 8-PSK(Phase Shift Keying), N/4 rows and 4 columns when the higher ordermodulation is 16-APSK (Amplitude Phase Shift Keying), and N/5 rows and 5columns when the higher order modulation is 32-APSK.
 3. A methodaccording to claim 1, wherein the encoder in the receiving stepgenerates the codeword according to a Low Density Parity Check (LDPC)code.
 4. A method according to claim 3, wherein the parity check matrixof the LDPC code is structured by restricting a triangular portion ofthe parity check matrix to zero values.
 5. A method according to claim3, wherein the higher order constellation represents a Quadrature PhaseShift Keying (QPSK) modulation scheme, the method further comprising:determining an i^(th) QPSK symbol based on the set of 2 i ^(th) and (2i+1)^(th) LDPC encoded bits, wherein i=0, 1, 2 . . . , N/2−1, and N isthe coded LDPC block size.
 6. A method according to claim 3, wherein thehigher order constellation represents an 8-PSK modulation scheme, themethod further comprising: determining an i^(th) 8-PSK symbol based onthe set of (N/3+i)^(th), (₂N/3+i)^(th) and i^(th) LDPC encoded bits,wherein i=0, 1, 2, . . . , N/3−1, and N is the coded LDPC block size. 7.A method according to claim 3, wherein the higher order constellationrepresents a 16-APSK (Amplitude Phase Shift Keying) modulation scheme,the method further comprising: determining an i^(th) 16-APSK symbolbased on the set of (N/2+2i+1)^(th), 2 i ^(th), (N/2+2i+1)^(th) and (2i+1)^(th) LDPC encoded bits, wherein i=0, 1, 2, . . . , N/3−1, and N isthe coded LDPC block size.
 8. A method according to claim 3, wherein thehigher order constellation represents a 32-APSK (Amplitude Phase ShiftKeying) modulation scheme, the method further comprising: determining ani^(th) 32-APSK symbol based on the set of (N/5+i)^(th), (2N/5+i)^(th),(4N/₅+i)^(th), (3N/5+i)^(th) and i^(th) LDPC encoded bits, wherein i=0,1, 2, . . . , N/5−1, and N is the coded LDPC block size.
 9. Acomputer-readable medium bearing instructions for transmitting encodedsignals, said instruction, being arranged, upon execution, to cause oneor more processors to perform the method of claim
 1. 10. A transmitterfor generating encoded signals, the transmitter comprising: an encoderconfigured to transform an input message into a codeword represented bya plurality of set of bits; and logic configured to map non-sequentiallyone set of bits into a higher order constellation, wherein a symbol ofthe higher order constellation corresponding to the one set of bits isoutput based on the mapping.
 11. A transmitter according to claim 10,wherein the N encoded bits are written to a block interleaver column bycolumn and read out row by row, and the block interleaver has N/3 rowsand 3 columns when the higher order modulation is 8-PSK (Phase ShiftKeying), N/4 rows and 4 columns when the higher order modulation is16-APSK (Amplitude Phase Shift Keying), and N/5 rows and 5 columns whenthe higher order modulation is 32-APSK.
 12. A transmitter according toclaim 11, wherein the encoder generates the codeword according to a LowDensity Parity Check (LDPC) code.
 13. A transmitter according to claim12, wherein the parity check matrix of the LDPC code is structured byrestricting a triangular portion of the parity check matrix to zerovalues.
 14. A transmitter according to claim 12, wherein the higherorder constellation represents a Quadrature Phase Shift Keying (QPSK)modulation scheme, and the logic is further configured to determine ani^(th) QPSK symbol based on the set of 2 i ^(th) and (2 i+1)^(th) LDPCencoded bits, wherein i=0, 1, 2 . . . , N/2−1, and N is the coded LDPCblock size.
 15. A transmitter according to claim 12, wherein the higherorder constellation represents an 8-PSK modulation scheme, and the logicis further configured to determine an i^(th) 8-PSK symbol based on theset of (N/3+i)^(th), (2N/3+i)^(th) and i^(th) LDPC encoded bits, whereini=0, 1, 2, . . . , N/3−1, and N is the coded LDPC block size.
 16. Atransmitter according to claim 12, wherein the higher orderconstellation represents a 16-APSK (Amplitude Phase Shift Keying)modulation scheme, and the logic is further configured to determine ani^(th) 16-APSK symbol based on the set of bits (N/2+2i)^(th), 2 i ^(th),(N/2+2i+1)^(th) and (2 i+1)^(th) LDPC encoded bits, wherein i=0, 1, 2, .. . , N/3−1, and N is the coded LDPC block size.
 17. A transmitteraccording to claim 12, wherein the higher order constellation representsa 32-APSK (Amplitude Phase Shift Keying) modulation scheme, and thelogic is further configured to determine an i^(th) 32-APSK symbol basedon the set of bits (N/5+i)^(th), (2N/5+i)^(th), (4N/5+i)^(th),(3N/5+i)^(th) and i^(th) LDPC encoded bits, wherein i=0, 1, 2, . . . ,N/5−1, and N is the coded LDPC block size.
 18. A method for processingencoded signals, the method comprising: demodulating a received encodedsignal representing a codeword, wherein the encoded signal beingmodulated according to a non-sequential mapping of a plurality of bitscorresponding to the codeword; and decoding the codeword associated withthe encoded signal.
 19. A method according to claim 18, wherein the Nencoded bits are written to a block interleaver column by column andread out row by row, and the block interleaver has N/3 rows and 3columns when the higher order modulation is 8-PSK (Phase Shift Keying),N/4 rows and 4 columns when the higher order modulation is 16-APSK(Amplitude Phase Shift Keying), and N/5 rows and 5 columns when thehigher order modulation is 32-APSK.
 20. A method according to claim 19,wherein the decoding step is according to a Low Density Parity Check(LDPC) code.
 21. A method according to claim 20, wherein the paritycheck matrix of the LDPC code is structured by restricting a triangularportion of the parity check matrix to zero values.
 22. A methodaccording to claim 20, wherein the higher order constellation representsa Quadrature Phase Shift Keying (QPSK) modulation scheme, and an i^(th)QPSK symbol is determined based on the set of 2 i ^(th) and (2 i+1)^(th)LDPC encoded bits, wherein i=0, 1, 2 . . . , N/2−1, and N is the codedLDPC block size.
 23. A method according to claim 20, wherein the higherorder constellation represents an 8-PSK modulation scheme, and an i^(th)8-PSK symbol is determined based on the set of (N/3+i)^(th),(2N/3+i)^(th) and i^(th) LDPC encoded bits, wherein i=0, 1, 2, . . . ,N/3−1, and N is the coded LDPC block size.
 24. A method according toclaim 20, wherein the higher order constellation represents a 16-APSK(Amplitude Phase Shift Keying) modulation scheme, and an i^(th) 16-APSKsymbol is determined based on the set of bits (N/2+2i)^(th), 2 i ^(th),(N/2+2i+1)^(th) and (2 i+1)^(th), LDPC encoded bits, wherein i=0, 1, 2,. . . , N/3−1, and N is the coded LDPC block size.
 25. A methodaccording to claim 20, wherein the higher order constellation representsa 32-APSK (Amplitude Phase Shift Keying) modulation scheme, and ani^(th) 32-APSK symbol is determined based on the set of bits(N/5+i)^(th), (2N/5+i)^(th), (4N/5+i)^(th), (3N/5 encoded bits, whereini=0, 1, 2, . . . , N/5−1, and N is the coded LDPC block size.
 26. Acomputer-readable medium bearing instructions processing encodedsignals, said instruction, being arranged, upon execution, to cause oneor more processors to perform the method of claim 18.